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Two replicas of spatially extended chaotic systems synchronize to a common spatio-temporal 
chaotic state when coupled above a critical strength. As a prototype of each single spatio-temporal 
chaotic system a lattice of maps interacting via power-law coupling is considered. The synchro- 
nization transition is studied as a non-equilibrium phase transition, and its critical properties are 
analyzed at varying the spatial interaction range as well as the nonlinearity of the dynamical units 
composing each system. In particular, continuous and discontinuous local maps are considered. In 
both cases the transitions are of the second order with critical indexes varying with the exponent 
characterizing the interaction range. For discontinuous maps it is numerically shown that the tran- 
sition belongs to the anomalous directed percolation (ADP) family of universality classes, previously 
identified for Levy-flight spreading of epidemic processes. For continuous maps, the critical expo- 
nents are different from those characterizing ADP, but apart from the nearest-neighbor case, the 
identification of the corresponding universality classes remains an open problem. Finally, to test 
the influence of deterministic correlations for the studied synchronization transitions, the chaotic 
dynamical evolutions are substituted by suitable stochastic models. In this framework and for the 
discontinuous case, it is possible to derive an effective Langevin description that corresponds to that 
proposed for ADP. 

PACS numbers: 05.45.Xt,05.45.-a,47.27.-i 



Synchronization is ubiquitous in Nature: neu- 
ronal populations, cardiac pacemakers, Joseph- 
son circuits and even coupled chaotic systems can 
synchronize during their activity. Remarkably, 
all these different phenomena can be described 
within a common framework represented by non- 
linear dynamics 

Synchronization of spatially extended chaotic 
systems is particularly interesting as it allows 
for transferring concepts and methods borrowed 
from statistical mechanics to nonlinear dynam- 
ics. Indeed, a direct connection between the crit- 
ical properties of non-equilibrium phase transi- 
tions and synchronization processes has been es- 
tablished [2, 3]. Noticeably, for diffusively cou- 
pled chaotic systems only two universality classes 

Directed Percolation and Multiplicative Noise 

encompass all the synchronization transitions. 
Originally these classes have been identified in 
completely different contexts such as epidemics 
spreading and pinning/depinning of interfaces 
to/from a substrate P, o]. The parallel between 
non-equilibrium critical phenomena and synchro- 
nization processes is possible only thanks to the 
erratic nature of the synchronized state, where 



chaos somehow mimics the presence of thermal 
noise in real systems. 

In this paper we propose an extension of such 
an analogy to synchronization transitions of long- 
range coupled systems. Long range interactions 
naturally appears in many circumstances rang- 
ing from neuronal networks to solid state physics. 
Therefore, a characterization of these transitions 
is definitely worth and of interest for a wide sci- 
entific community. 

I. INTRODUCTION 

Since its discovery chaotic synchronization is a 
very active and important field of research |1|, with 
applications in such diverse fields as secure communi- 
cations_0, [1], semiconductor lasers @, chemical reac- 
tions living systems [Tl| and electrically coupled 
neurons in vitro |12l |. 

The phenomenology of the synchronization transition 
(ST) is particularly rich for spatially extended chaotic 
systems, where a parallel with critical phenomena can 
be drawn [l3|. In the last decade, an ongoing research 
activity has been devoted to relate chaotic STs to non- 
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equilibrium phase transitions 

Ho, HH m, [2J, Hi]. Nowadays, it is well established that 
the synchronization of two replicas of a chaotic system 
coupled via spatio-temporal noise or via local trans- 
verse coupling j3|, can be characterized as a continuous 
non-equilibrium transition from an active to an absorbing 
phase. In systems with nearest-neighbor (NN) couplings, 
it has been found that the STs belong to two universal- 
ity classes, depending on the local dynamical features, 
namely: directed percolation (DP) [1, [2^ or multiplica- 
tive noise (MN) These studies have been mainly 
performed for coupled map lattices (CML), which are 
prototype models for systems exhibiting spatio-temporal 
chaos [2^1 . An order parameter for the STs is represented 
by the propagation velocity of information Vi [28L [2^ : in 
the synchronized (desynchronized) state this quantity is 
zero (finite). Moreover, for continuous maps, where the 
dynamics is dominated by linear effects Vi vanishes to- 
gether with the transverse Lyapunov exponent. In this 
case the critical indexes associated to the ST are of the 
MN kind On the other hand, for discontinuous (or 
quasi-discontinuous [131 ) maps, due to the predominance 
of nonlinear effects, the replicas synchronize for defi- 
nitely negative transverse Lyapunov exponents 0, [29| . 
For these maps the ST is now characterized by critical 
properties typical of DP, analogously to what happens 
for cellular automata [l^, HQi] . These findings have been 
also confirmed by the analysis of stochastic models mim- 
icking the synchronization of continuous/discontinuous 
maps 
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Recently, synchronization has been also studied for 
chaotic systems presenting long-range interactions, which 
are relevant to man y re al contexts, such as disease spread 
via aviation traffic |3l| , neuron populations [s^ l, J oseph- 
son junctions [s^ and cardiac pacemaker cells [3^ . 

Typically, long-range interactions are introduced by 
considering CMLs with coupling decaying as a power- 
law [U, [3^, m, [13] ■ In such a model, when the local 
dynamics is dominated by nonlinear mechanisms, the 
STs have been shown to belong to a family of univer- 
sality classes known as Anomalous Directed Percolation 
(ADP) '24, 38]. ADP has been previously identified for 
epidemic spreading whenever the infective agent can per- 
form unrestricted Levy flights Jsl, Such processes, 
originally introduced in Ref. [4l|, can be modeled by 
assuming, e.g. in d = 1, that the disease propagates 
from an infected site to any other with a probability 
P{r) ^ 7-~(^+'^' algebraically decaying with the spatial 
distance r, where a controls the interaction range. Hin- 
richsen and Howard (HH) have numerically shown for a 
stochastic lattice model (generalizing directed bond per- 
colation ) th at the critical exponents vary continuously 
with a [40|. These findings confirm previous theoreti- 
cal results [11] indicating that usual DP should be re- 
covered for sufficiently short-ranged coupling (namely, 
(T > (Tc = 2.0677(2)) and that a mean-field description 
should become exact for a < am = 0.5 (for a recent and 
exhaustive review see Ref. 1381). 



In this paper, we first reconsider, by performing more 
accurate estimations of the critical exponents, the results 
obtained for NN interactions both for deterministic as 
well as for stochastic local dynamics [l^ [13] . Then we 
focus on the synchronization transition in systems with 
power-law coupling. In this context we extend previ- 
ous analysis for discontinuous deterministic maps [13] 
to a stochastic version of the model, and we confirm 
that whenever the STs are driven by nonlinear effects 
the associated critical properties belong to the anoma- 
lous directed percolation universality classes. Finally, we 
present a first evaluation of the critical properties of the 
synchronization transitions for continuous chaotic maps 
with power-law interactions. 

The paper is organized as follows: the employed de- 
terministic and stochastic models are introduced in the 
next section, while Section HITl is devoted to the methods 
used to study their critical properties. In Section lTVl STs 
are re-examined for nearest-neighbor interacting systems; 
while Sec. |V] is focused on power-law coupled systems. 
Conclusions and perspectives are discussed in Sec. I VII 
The Appendix presents a derivation of the field equa- 
tion, known to reproduce ADP critical phenomena, for 
the stochastic model with power-law interactions. 



II. MODELS 

In this paper we investigate the synchronization be- 
tween two replicas of CMLs transversally coupled accord- 
ing to the following scheme 



xdt + 1) = {l-j)F{£,{t))+jF{y,{t)) 
y,it + l) = (l-7)F(y~(t))+7F(f,(t)) 



(1) 



where i — 1, . . . ,L is the discrete spatial index with L 
denoting the system size; Xi (t) , yi (t) G JO : 1] are the state 
variables and F{x) determines the local dynamics on each 
site of the lattice (see below for its specification). The 
parameter 7 sets the strength of the site-by-site coupling 
between the two replicas. 

The variables Zi € {xi,yi} indicate spatially averaged 
quantities. In particular, two different kind of spatial av- 
erages are introduced to reproduce short and long-ranged 
interactions. The former represents a discretized version 
of spatial diffusion among nearest neighbour sites and 
reads: 



■Si = (1 - f-)zt + 7r(.2i-i + 2j+i) , 



(2) 



where e measures the intensity of the diffusive interaction 
within each replica. The latter is obtained b y co nsidering 
a coupling decaying as a power law [U, [IE 113] > i-e. 



L' . 

Zi = e)2j -I- 37^ 2^ 



fc=i 



(3) 



where a tunes the interaction range (the rationale for 
defining the exponent as 1 + a will become clear in 
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Sec. |V]). For cr oo, Eq. ([3]) reduces to the diffusive 
case while for cr = —1 it corresponds to the mean- 
field coupling, usually employed in the study of glo bally 
coupled maps [43|. Actually, as discussed in Ref. [2l| . 
the large scale properties of the system should coincide 
with those at ct = for any a E [—1:0]. Since the sum 
extends up to L' = (L — l)/2 the model is well defined 

only for odd L- values, and ri{a) — 2 ^^^j^ fc~(^+'^^ is a 
normalization factor to ensure that z g [0 : 1]. For both 
kinds of coupling ^ and ^ periodic boundary condi- 
tions are assumed and the diffusive coupling e is set to 
2/3. Notice that the results do not seem to depend on 
the chosen value of e. 

Several studies 0, S, [13, HH HI have shown that the 
system dynamics is strongly influenced by the continu- 
ity properties of the local map F{x) [23|. In particu- 
lar, the spatio-temporal propagation of disturbances [111 
and the properties of synchronization [^ are noticeably 
different if the map is continuous or (quasi) discontin- 
uous. For this reason, throughout this work, we con- 
sider two different deterministic maps: the Bernoulli map 
F(x) = 2x(modl), which is characterized by a discon- 
tinuity making dominant the nonlinear effects, and the 
tent map F{x) = 1 — 2\x — 1/2|, whose dynamical prop- 
erties are well captured by linear analysis. 

The synchronization transition is observed while in- 
creasing the coupling 7: above a critical value 7c the 
two replicas synchronize onto the same chaotic trajec- 
tory. In other words, for 7 > 7c the synchronization 
error Wi{t) = \xi{t) — yi{t)\ vanishes for sufficiently long 
times. 

In parallel with Eq. ([1]), we also consider a stochas- 
tic model which mimics the evolution of the synchro- 
nization error (difference field) Wi{t) in proximity of the 
ST. Stochastic dynamics is expected to deplete the role 
played by spatio-temporal correlations, unavoidable in 
deterministic systems, and thus should allow more accu- 
rate estimations of the critical exponents. 

In particular, we consider a stochastic model able to 
reproduce the main features of ST both for continu- 
ous and discontinuous maps. We consider the so-called 
Random Multiplier (RM) model, originally proposed in 
Refs. (B 117|. The quantity Wi{t) defined in each site 
of the chain evolves according to the following stochastic 
dynamics: 

^ ' [ awi[t), w.p. 1—p 

(4) 

^ ' [ awi[t), w.p. 1—p 

(5) 

where Wi (t) g [0 : 1] and Wi indicates the spatial averages 
that depending on the nature of the coupling, NN or long- 
range, is given by or The parameter A controls 
the nonlinear effects, while a the linear ones. For the local 
maps considered in this article a can be identified with 



2(1 — 27) [T3|. Therefore, to vary the coupling 7 in the 
deterministic model amounts to modify a in the stochas- 
tic one. Notice that, due to the relationship between the 
parameters 7 and a, in the RM model to achieve syn- 
chronization a should be reduced. For sufRciently small 
(resp. large) A the model mimics the dynamics 

of (quasi)-discontinuous (resp. continuous) maps, where 
nonlinear (resp. linear) effects are predominant. Follow- 
ing Ref. [141 , in order to reproduce the behaviors of the 
Bernoulli or of continuous-like maps we fixed A = or 
A = 0.2, respectively. 



III. CHARACTERIZATION OF THE 
SYNCHRONIZATION TRANSITIONS 

The phenomenology of the synchronization transitions 
at varying the coupling range and the type of local dy- 
namics can be visualized by looking at Fig. [1] The figure 
displays the difference field Wi{t) for various values of 
a with transversal coupling just above the synchroniza- 
tion transition (7 > 7c) for Bernoulli as well as for tent 
maps. The spatio-temporal evolution of Wilt) is strongly 
dependent not only on the continuous / discontinuous na- 
ture of the local dynamics, but also on the interactions 
range within each replica. For the discontinuous maps 
(Fig. [T^-d), percolating structures are clearly visible in 
the short range limit (i.e. large a). However, as the 
range of the interaction increases (i.e. a decreases) these 
spatial structures tend to be smoothed out and, finally, 
for cr < 0.5, they are no more detectable. When the local 
dynamics is governed by a continuous map, the qualita- 
tive results change completely (FiglT^-h). In this case, 
desynchronization (resurgence) phenomena are possible 
within the already synchronized areas and this leads to a 
less defined distinction between synchronized and desyn- 
chronized regions. However, a disappearance of the spa- 
tial structures is once more observable by increasing the 
interaction range. 

The patterns observed in Fig. [T] suggest a link be- 
tween synchronization transitions and non-equilibrium 
phase transitions from an active to an absorbing (quasi- 
absorbing) phase [1, [1|. In particular. Fig. [T^ recalls the 
patterns observed for Directed Percolation (DP) [^. [25j. 
a contact process usually employed in the description of 
epidemics spreading. DP has been the subject of active 
research in the last 20 years and, notwithstanding ex- 
act analytical results are still lacking, detailed numerical 
studies have been able to determine its critical proper- 
ties up to a noticeable accuracy [i^ . Quite recently, also 
experimental measurements of the critical DP indexes 
have been reported for a transition to spatio-temporal 
intermittency in a quasi one-dimensional system consist- 
ing of a ring of ferrofluidic spikes 0| , and for a transi- 
tion between two topologically different turbulent states 
in a quasi two-dimensional layer of nematic liquid crys- 
tals [45| . Fig. [TJj resembles instead surface roughening, 
indeed Kurths and Pikovsky A6] have proved that within 





FIG. 1: (Color online) The spatio-temporal evolution of the synchronization error for power-law coupled maps is depicted for 
local dynamics given by the Bernoulli map (upper row) and the tent map (lower row). Each column corresponds to different 
values of the exponent a: namely, ct = 5 (panels (a) and (e)); cr = 2 (panels (b) and (f)); a = 1.2 (panels (c) and (g)) and 
o — 0.5 (panels (d) and (h)). The horizontal axis represents space, while the vertical one, time. Note that the patterns observed 
for NN interaction closely resembles those reported in (a) and (e) for discontinuous and continuous maps, respectively. The 
color code represents the synchronization error (in logarithmic scale), ranging from white {wi{t) —> 0) to black {wi{t) — 1), 
larger values correspond to darker tones. The results reported in the panels have been obtained by employing the modified 
model (flit for g = 4 and M = 5. 



a linear framework the evolution of the difference field 
of coupled CMLs is described by a Kardar-Parisi-Zhang 
equation [i^] in the presence of a hard wall. This lat- 
ter equation can be put in directed relationship with 
the Multiplicative Noise (MN) Langevin equation Q, 
which reproduce critical properties of non-equilibrium 
pinning^lcmnning as well as wetting-dewctting transi- 
tions [ig. l-jsj. 

As mentioned in the Introduction, the above analo- 
gies are not merely qualitative, indeed several stud- 
ies have quantitatively shown that, depend- 
ing on the prevalence of linear (resp. nonlinear) effects 
in the local dynamics, ST in diffusively coupled systems 
belongs to the MN (resp. DP) universality class. In par- 
ticular, by changing the value of the local multiplier (a 
control parameter), one can pass continuously from one 
class to the other flT'], thus suggesting that both these 
non-equilibrium transitions can be described within a sin- 
gle field-theoretic framework. This is confirmed by fur- 
ther results reported in Ref. [l^ for a microscopic model 
of wetting transitions and in Ref. [l^ for KPZ equation 
with an attractive wall. Indeed, for such processes the 
time evolution of the local density of active sites n(x,t) 
admits a formally quite similar Langevin description : 



dtn = DnV^u + Tn~ TnP + g{n)S, . 



(6) 



Such equation recovers the Reggeon-Field Theory (RFT) 
describing DP for p — 2 and g{n) oc ^/n, while the mini- 
mal model for MN can be obtained from Eq. ([6]) whenever 
g{n) oc n. In the synchronization context, n represents 
the coarse grained synchronization error obtained by av- 
eraging Wi{t) = \xi{t) — yi{t)\ over a suitable space-time 



cell. The parameters entering Eq. ^ are the diffusion 
coefficient Dn (corresponding to e); r that measures the 
distance from the critical point (i.e. from the synchro- 
nization threshold 7c) and the amplitude of the nonlinear 
term T. Finally, ^ is a zero-average ^-correlated (in space 
and time) Gaussian noise field with unit variance. 

In Ref. 39, 40] the RFT has been extended to include 
long-range interactions, leading to the following stochas- 
tic equation 

dtn = DW'^n + Da^^u + tu - Fn^ -t- g{n)£, , (7) 

which generalizes ([6]) through the addition of an anoma- 
lous diffusion term with coefficient Da and range of in- 
teraction parametrized by a. For g{n) oc y^, Eq. ([7|) 
describes the ADP universality class while, to the best of 
our knowledge, such an equation has never been studied 
in the context of anomalous MN; i.e for g(n) oc n. 

Recently in Ref. [l^l we have numerically shown that 
STs observed for model ^ with power law coupling 
([5]) and equipped with discontinuous maps (namely, 
Bernoulli maps) belong to ADP universality class. In 
this context the extension of the analysis to the stochas- 
tic models (j4l5p would be particularly interesting, for two 
reasons. Firstly, the model (I4I5P is known to reproduce 
quite well the phenomenology of synchronization for sys- 
tems like ([T]) with short-range coupling [13]. Second, after 
a suitable coarse-graining, such models for A = are ef- 
fectively described by the Langevin equation derived for 
RFT, providing further support to the connection with 
DP [13]. Here, in the Appendix, we show that this kind 
of mapping can be extended also to the power-law in- 
teracting model and, in this case, Eq. ([7]) is recovered. 
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Furthermore, the aim of the present work is to extend 
the studies reported in Ref. [2J] also to continuous maps. 

Before reporting a detailed analysis of the synchroniza- 
tion transition for the above defined models, let us in- 
troduce the exponents employed to characterize STs as 
non-equilibrium transitions. These are defined through 
the spatially averaged synchronization error p-y{t) = 
'^iWi{t)/ L which, at sufficiently long times, vanishes 
whenever a complete synchronization is achieved, i.e., 
p* = limt^oo Pjit) = for 7 > whilst it remains 
finite at any time in the desynchronized state, i.e. for 
7 < 7c- 

The order parameter py (t) allows us to define the crit- 
ical exponents 9, (3 and z which characterize the transi- 
tion: 9 is the exponent that rules the temporal scaling of 
the order parameter at the critical point 7 = 7c, namely 

P,At)-t''- (8) 

For 7 < 7c, /O* 7^ and one has 

p; ~ (7c - if (9) 

which defines the critical exponent (3. Finally, the dy- 
namical exponent z can be defined in terms of the finite- 
size scaling relation, valid at the critical point: 

PySt)-L'''f{tlL'). (fO) 

These three indexes are sufficient to fully characterize the 
transition, since from their knowledge all the other criti- 
cal exponents can be derived. In particular, the exponent 
(resp. vw) ruling the divergence of spatial (resp. tem- 
poral) correlation lenght at the critical point is given by 

= ^ (resp. = ^ ) (ff) 

A detailed description of the numerical estimation of 9, 
(3 and z is reported in the following sections. 

IV. SHORT-RANGE INTERACTIONS 

Let us now consider the synchronization of two replicas 
of diffusively coupled map lattices. We start with the 
case of coupled tent maps where the linear analysis is 
sufficient to identify the critical point and we can thus 
limit to study the behavior of the difference field 'Wi{t) 
in the tangent space. In particular, the norm of this 
field grows at a rate given by the transverse Lyapunov 
exponent [3]: 

= lim =ln(I-27)+A; (12) 

t — ^00 

where A is the maximal Lyapunov exponent of a single 
CML. The synchronization transition occurs at the point 
where Aj^ vanishes, locating the critical coupling to 

7c = ^(l-e-^). (13) 



For this model, it has been suggested that 7c = 
0.17615(5) [3], by extrapolating from finite size mea- 
surements the asymptotic value for A. We have es- 
timated the synchronization value from the scaling of 
the density p-y{t) for chains of length L = 2^^ finding 
7c — 0.17616(2), in agreement with the previous estima- 
tion. The corresponding critical exponents are reported 
in Table HI we observe that these values are in reasonable 
agreement with the one reported in literature for the MN 
universality class [1, E^, apart for 9, which is larger than 
the corresponding MN value. This difference is proba- 
bly due to the presence of correlations in the dynamical 
evolution of these maps, that limits the time range over 
which a critical scaling of the density ^^(i) can be found. 
In particular, in our analysis as well as in Ref. |3*|, the 
time interval where it is possible to observe a clear power- 
law scaling is limited to [10^ : 4.5 x 10^], since at longer 
times the density (above and below 7c) always saturates 
to a (small) constant value. As we show in the follow- 
ing, much more accurate estimations can be obtained by 
replacing the tent map with its stochastic version, thus 
reinforcing the hypothesis that scaling laws are hindered 
by long-time dynamical correlations. 

For coupled Bernoulli maps linear analysis fails in lo- 
cating the synchronization transition, therefore we are 
forced to directly investigate the scaling of the density 
Pj{t). In particular, for L = 2^"^ we have estimated a 
critical value 7c = 0.28752(1), in agreement with Ref. 3 
(see Table U for the values of the critical exponents). In 
this case finite size effects and time correlations are less 
relevant, since good scaling can be already observed in 
the time interval [10^ : 10^] for much shorter chains. 

Let us now turn our attention to the Random Multi- 
plier model (I4l5p . which reproduces the Bernoulli map 
behavior for A = and that of a generic continuous map 
for A = 0.2 [ll,[i3l. 

For the RM model with A = the scaling laws are 
reported in Fig. [2^-c. In this case the critical point is at 
ac = 0.60615(5) and the measured indexes are in perfect 
agreement with the best estimates of DP's exponents. 
9 and (3 coincide with the DP values up to the fourth 
and third digit, respectively. The evaluation of z is less 








P 


z 


Bernoulli 


0.159(1) 


0.27(1) 


1.58(4) 


RM (A = 0) 


0.1595(4) 


0.276(2) 


1.58(2) 


DP [43] 


0.159464(6) 


0.276486(6) 


1.580745(6) 


Tent 


1.275(15) 


1.70(8) 


1.5(1) 


RM (A = 0.2) 


1.13(3) 


1.67(3) 


1.53(6) 


MN [49] 


1.10(5) 


1.70(5) 


1.53(7) 


MN [50] 


1.184(10) 


1.776(15) 





TABLE I: Critical exponents for nearest-neighbor coupled 
tent and Bernoulli maps as well as for the corresponding 
stochastic models. Also the best estimations of the critical 
indexes for DP and MN universality classes are reported. 
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FIG. 2: (top) Scaling relationships for the critical exponents for the RM model (|4l5p with NN coupling for A = 0: (a) pa{t) vs 
t is displayed for a — 0.60620, 0.60615, 0.60610 with L — 2^® maps and averaging over 5 x 10^ different realizations, the inset 
shows the logarithmic derivative dln{pa{t))/d\n{t); (b) various curves Pa(t) vs t are reported for a < Oc revealing the saturation 
to p*; in the inset p* vs (a — Uc) is reported together with the best fit giving the /3 estimation, here L = 2^® maps and the 
average is over 500 realizations; (c) data collapse based on (|10|) . averaging over lO'^-lO^ initial conditions, (bottom) Same as 
top for the RM model with A = 0.2: (d) analogous to (a) for a = 0.566938, 0.566942, 0.566945 with L = 5 x lO'^ and averaging 
over 1 — 2 X 10^ realizations, here in the inset we report pac(t)t^ vs t; (e) same as (b) with L = 5 x 10^ and averaging over 
2^ — 1.2 X 10^ initial conditions; (f) analogous to (c), averaging over 10^ initial conditions. The values of the estimated critical 
exponents are reported in Table |I] 



accurate, since it relies on a data collapse of finite size 
estimates (|10p . 

In the case A — 0.2, by considering the scaling in time 
of the density ([8]) for L = 5 x lO'^ maps over a time 
span t € [lO'^ : 10^], we have found that the critical 
point is located at Qc = 0.566942(4). In particular, 9 
has been estimated by a best fit to a power-law over 
more than two decades (Fig. [2]i) . The quality of the es- 
timation of the other two exponents can be appreciated 
from Fig. [5^ and Fig. [2f. The values of the critical in- 
dexes are in fairly good agreement with those reported 
in Ref. [49[, obtained by considering a time and space 
discretized version of the corresponding Langevin equa- 
tion, i.e. Eq. ^ with g{n) cx ^/n. It should be remarked 
that the exponents recently reported in Ref. [s^, for a 
suitable lattice model reproducing KPZ-type interface 
growth, are slightly larger than ours, while the 0- value 
(namely, 9 = 7/6 = 1.1666) conjectured by Droz and 
Lipowski 20'] is, within the error bars, consistent with 
our estimation. 



V. LONG-RANGE INTERACTIONS 

In the presence of long-range interactions, accurate 
numerical analysis of the synchronization transition re- 



quires huge time costs, this can be appreciated by notic- 
ing that at each time step 0{L^) operations are involved 
(see the coupling definition in ([3])). Simulation times be- 
come prohibitive as the range of the interaction increases 
(i.e. the exponent cr decreases), because finite size effects 
becomes more relevant and larger sizes are needed. It is 
therefore fundamental to reduce the CPU cost for each 
time step. As in Refs. [l^, H^, to achieve such aim we 
consider the following coupling scheme ([3]): 



a-e)z, 



ri{a) 



(Zi_l + Zi+l) + 



M 

E 



m=l 



ijm{q)Y 



(14) 



where jm{q) 



1 and M = log„(L/2) |5li; the 



model with full coupling ([3]) is recovered for jmil) = rn, 
M — L — 1 and by substituting a — > 1 + a (see below for 
more details) . Clearly, the new choice is very convenient 
because each updating step only requires 0{L log^ L) op- 
erations instead of 0{L?) needed for the fully coupled 
case. The parameter q is typically chosen as q — 2, 4 
and 8, while the normalization factor is given by 77(17) = 

2[l + E™=l,M0m(9));"]"'- 

The critical properties of the reduced model {jm{q) = 
q™ — 1) with exponent a map into those of the fully cou- 
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FIG. 3; (Color online) Critical exponents (a) 6, (b) (3 and (c) z for coupled Bernoulli maps and RM models (for A = 0) as a 
function of the long-range power exponent cr compared with those obtained by HH in Ref. [i^l for a stochastic lattice model. 
Symbols are explained in the legend. The dotted horizontal line represents the best estimate of the critical exponents for usual 
DP, i.e. for nearest-neighbor spreading. 



HH 



B g = 2 



B g = 4 



B q = 8 



RM g = 4 



/3 



f3 z 



0.2 
0.5 
0.6 
0.8 
1.0 
1.1 
1.2 
1.4 
1.5 
1.6 
1.8 
2.0 
2.2 
2.4 
2.5 
3.0 
5.0 



0.99(4 
0.94(4 
0.86(4 
0.67(3 
0.52(3 

0.40(3 
0.33(3 

0.27(3 
0.24(4 
0.21(4 
0.19(4 
0.17(4 



0.99(4 
0.95(6 
0.88(5 
0.76(4; 
0.65(3 

0.56(3 
0.49(3 

0.43(3 
0.39(3 
0.34(3 
0.32(3 
0.30(3 



0.21(2 
0.54(2 
0.58(2 
0.71(2 
0.82(2 

0.96(3 
1.09(3 

1.21(3 
1.32(3 
1.43(3 
1.49(4 
1.53(5 



0.50(1) 0.67(5) 0.83(4) 

0.35(1) 0.52(5) 0.96(4) 

0.30(1) 0.49(5) 1.14(4) 

0.19(2) 0.35(5) 1.40(3) 

0.16(1) 0.29(5) 1.58(4) 



0.73(2) 0.91(5) 

0.48(2) 0.64(5) 

0.40(5) 0.57(4) 0.92(5) 

0.33(5) 0.49(5) 1.05(5) 



0.20(2) 0.35(3) 1.45(5) 



0.97(4) 1.02(5) 0.50(5) 
0.68(5) 0.80(5) 0.67(5) 



0.52(1) 0.65(1) 0.75(5) 

0.33(1) 0.48(2) 1.06(4) 
0.27(1) 0.41(1) 1.20(4) 



0.159538 0.276 1.58 



TABLE II: Summary of the results for the critical exponents (also shown in Fig. |3]) corresponding to coupled Bernoulli maps 
and to chains of RM-models (with A = 0) at various values of a. Results taken from Ref. [4(| by HH are also reported for 
comparison. 



pled one (j™ —m) with exponent c/c = cr+l. This can 
be easily understood by noticing that the two versions 
of the model should display the same critical behav- 
ior once the spatial interactions scale analogously. For 
the modified model the coupling weight over the inter- 
val \jm{fl) ■ jm+iiq)], containing a single coupled site, is 
simply given by 



1 



1 



Jm+l 



(m+l)(T 



while over the same interval the weight for the fully cou- 
pled model can be estimated by evaulating the sum 



1 



1 



fc=Jm(9) + l 



(m+l)(<T^,-l) 



By comparing the two expressions it is thus clear that the 
two weights scale in the same manner when afc = (J + l, 



24\ 



as we have also numerically verified in Ref. 

Before discussing the numerical results we observe that 
the coupling scheme (ITil) can be easily implemented both 
for coupled maps as well as for the stochastic model. 



A. Discontinuous maps and ADP universality class 

Synchronization transition of long-ranged coupled sys- 
tems with discontinuous maps has been previously stud- 
ied in Ref. [2l[, as far as the self-synchronization of a 
single chain is concerned, and in Ref. [24j, where the 
synchronization of two (transversely coupled) replicas is 
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discussed. Here we summarize the results obtained in 
Ref. lill for Bernoulli maps coupled with reduced cou- 
pling scheme [111 and extend the analysis to the stochas- 
tic RM model with A = in the same coupling condi- 
tions. 

In order to observe clean scaling laws, we have em- 
ployed the reduced scheme with increasingly larger q- 
values for smaller a. This amounts to chain lengths vary- 
ing from L = for a = 5.0 up to L = 2^^ for a = 0.5. 
Analogously to the NN case, due to the failure of linear 
analysis, the critical point 7c has been located by consid- 
ering the scaling relation ([5]), in particular from the crit- 
ical power-law decay the exponent 9 can be estimated. 
The critical index /3 can be derived by considering the 
asymptotic values p* below the transition (7 < 7c) as in 
([5]), finally z is estimated by data collapse based on the 
finite-size relation (fTO| at 7 = 7c. As shown in Ref. [23 |. 
the quality of the scaling behavior is very good. 

Table mi and Fig. [3] summarize our findings and present 
the new critical exponents estimated for the RM model 
(|11)-(IS|), where we have employed the coupling scheme 
(fTi)) with q = 4 and M = 7 (for <t = 1.4 and 1.6) and 
M = 9 (for a = 1.0). As one can appreciate from Ta- 
ble [IT] the estimated errors are smaller for the RM than 
for the Bernoulli map. In tabic HIl we also report the crit- 
ical indexes found by Hinrichsen and Howard for the 
Anomalous DP universality classes. As one can see the 
results obtained for the Bernoulli map approach those of 
HH for increasing q, i.e. for chains with larger sizes. This 
tendency suggests that the observed discrepancies are 
mainly due to finite size effects, which however seem to 
be less severe for the stochastic model. Indeed by consid- 
ering the RM model with chain sizes L ^ 3 x 10'^-5x 10^ 
we found that the exponents coincides, within the error 
bars, with those reported by HH. 

According to the analysis reported in Ref. [s^, H^j , the 
upper critical dimension dc is 2a and therefore the mean- 
field regime in one dimensional chain should establish be- 
low am = 0.5. The corresponding critical exponents are 
9m F = Pmf = 1 and zmf = f (in accordance the tem- 
poral and spatial correlation length exponents are given 
by i^Wj^jp — 1 and h'±MF — ^l<^\ Therefore, the long- 
range nature of the interactions is reflected in the mean- 
field regime only by the scaling of the spatial correlations, 
since for short-range interactions the exponents coincide 
apart zmf = 2 (and correspondingly z^_lm_f = 1/2). Our 
data reported at a = 0.5 for the Bernoulli maps confirm 
the mean-field expectations, but due to computational 
limitations we could not explore smaller a-valucs. 

For ADP, it has also been shown [s^ |4^ that the fol- 
lowing hyper-scaling relation holds 

d = l-a+{l-29)z = 0] (15) 

in the tr-range [cr„i ■ Cc]. The CTc value at which the be- 
haviour of the system should cross over to usual DP can 
be directly estimated by inserting the estimated values 
of the DP esponents in Eq. (fT5)l . Quite astonishingly the 
cross over takes place at CTc = 2.0677(2) > 2, as sug- 



gested also by field theoretic arguments [H, H^. As 
shown in Table Hill the relation is fulfilled within the 
error bars for the Bernoulli coupled maps as well as for 
its stochastic version. Notice that, 6 departs from zero 
only at a = 3.0 > CTc, where usual DP scalings are ex- 
pected. Finally, in the Appendix we show that, after a 
suitable spatio-temporal coarse-graining, the stochastic 
model can be effectively described by a Langevin equa- 
tion of the form ([7]) corresponding to the field theoretic 
description associated to ADP. These results further re- 
inforce the parallel between ST induced by nonlinear ef- 
fects, in the presence of power-law interactions, and ADP. 



B. Continuous maps 

Let us finally consider two coupled replicas with local 
dynamics given by the tent map. In this case the critical 
point can be estimated by the vanishing of the trans- 
verse Lyapunov exponent accordingly to the expression 
(1131) . Furthermore, we have independently evaluated 7c 





3.0 [2] 


2.0 [4] 


1.4 [4] 


1.2 [4] 


0.8 [8] 


0.5 [8] 


s 


0.93(6) 


-0.13(9) 


-0.04(12) 


-0.03(10) 


-0.04(9) 


0.03(9) 



TABLE III: Scaling relation (|15p for various values of a, for 
each measurement the corresponding basis q is reported. For 
the RM-model with q = 4 we obtained: for ct = 1.0 <5 = 
-0.03(2), for cr = 1.4 (5 = -0.04(4) and for cr = 1.6 5 = 
—0.05(4). The measured S are compatible with zero in the 
range [0.5 : 2.0677]. See text for details. 




0.5 1 L5 2 2.5 3 3.5 



a 

FIG. 4: 6 and as a function of a for the tent map with 
power-law coupling. The shaded upper (resp. lower) region 
corresponds to the best estimate of 6 (resp. /3) reported in 
literature [l^, [H^ for the MN class in the case of short 
range interactions. The thick upper (resp. lower) horizontal 
lines to the corresponding estimates for coupled tent maps 
with NN coupling reported in Table |T] The employed system 
sizes vary between L — 2^^ and L = 2^^ with q — A for the 
long-range models. For systems with NN interactions we used 
sizes L — 2^^ for the ^-estimation and L = 2^* for f3. 
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1.6 
1.4 
1.2 - 
1 




1 2 3 4 5 

o 

FIG. 5: z vs a for the tent map with power-law couphng. 
The shaded region corresponds the best estimate of z reported 
in Uterature [Sj for the MN class in the case of short range 
interactions. For the finite-size scaling chains of length L = 
2* — 2^^ have been employed. 



by examining the critical behaviour of the density pj (t) . 
Usually the Lyapunov approach was able to locate 7c 
within a precision of 1 — 2 x 10^* or even better. 

In the present case, to simplify the analysis we keep 
the parameter g = 4 and we analyze system sizes ranging 
from L = 2^^ to i = 2^^. The results for the exponents 
6 and /3 are shown in Fig. [4] while z in Fig. [5l we ob- 
serve that these critical indexes for ct > 3 tend to the 
results reported for the MN universality class (see Table 
|T]). For smaller cr- values and /3 (resp. z) increase (resp. 
decreases) and appear to saturate (resp. to vanish) for 
a < 1. While the numerical values are completely differ- 
ent from the ADP ones, the general trends are analogous 
to those observed for coupled Bernoulli maps. This sug- 
gests that in the present situation the mean-field value 
for the indexes are 9mf ~ 2.2(1), f3MF 2.6(1) and by 
assuming a linear dependence among z and the power- 
law exponent one finds that zmf ~ 0.65(5) x (cr — ctq), 
with (To ~ 0.40(5). At variance with the mean-field be- 
haviour of the ADP case, Omf and Pmf do not coincide, 
thus suggesting that t^\\j^p 1.18(10) and z seems to 
vanish already for a ~ 0.40(5). Finally, from the nu- 
merical estimates of z we can conjecture that i^±mf ^ 
1.8(5)/(cr — (To)- Apart from the poor quality of the esti- 
mates, it is important to stress that, analogously to ADP, 
the presence of long-range interactions reflects mainly in 
the scaling of the spatial correlations. The mean-field 
limit for systems with multiplicative noise has been ex- 
amined recently only for short-range couplings fs^ [ssj . 
The authors of the papers agree on the values of the cor- 
relation exponents i^Wj^.jp = 1 and v^mf — 1/2; which 
implies zmf = 2 and 6mf — Pmf- However, the value 
of 9mf seems to depend on the approximation employed 
to derive the mean-field. In particular, in H, m a 
value Omf = 5/3 is reported, while in [1, [5^ is predicted 



that 9mf is a nonuniversal scaling exponent dependent 
on noise amplitude, diffusion coefficient, and non-linear 
term. Therefore the mean-field analysis is still unclear 
for the multiplicative noise case with usual diffusion (as 
discussed also in ,5^] ) , moreover a study for the fractional 
diffusion case is still to be addressed. 

The analysis of the corresponding stochastic model, i.e. 
the RM model with A = 0.2 with power-law coupling, re- 
veal enormous fluctuations in the asymptotic behavior of 
the density induced by abrupt synchronizations of the 
whole chain. Thus making impossible any estimation of 
critical scaling laws, even by employing quite long chains 
(namely, L = 2^^). Our analysis cannot rule out the pos- 
sibility that for the RM model the transitions became 
discontinuous. Similarly to what done in Ref. (ssj . a 
detailed analysis of the scaling of the synchronized clus- 
ters with the system size should be performed to address 
this point, but this goes beyond our scopes. Moreover, 
it is quite astonishing that the analysis of the determin- 
istic case was somehow clearer. We can conjecture that 
the correlations induced by the deterministic evolution 
of coupled tent maps prevent the abrupt synchronization 
of large islands within the chain. However, this is just a 
working hypothesis to be investigated in the future. 



VI. CONCLUSIONS 

The critical properties of the synchronization transi- 
tions among replicas of chaotic and stochastic spatially 
extended systems have been numerically estimated both 
for diffusively coupled and for power-law interacting sys- 
tems. In particular, we focus on the differences between 
transitions dominated by linear and non-linear mecha- 
nisms. 

For nearest-neighbour coupling, our analysis confirm 
previous findings indicating that the transitions are al- 
ways continuous, while two distinct universality classes 
characterize the transition depending on the nature of 
the local dynamics. For continuous (resp. discontinuous) 
local maps the critical properties correspond to those of 
Multiplicative Noise (resp. Directed Percolation) non- 
equilibrium phase transitions. 

The introduction of a power-law coupling modifies the 
critical properties of the STs, in particular for (quasi- 
) discontinuous maps all the studied transitions can be 
gathered in a unique family of universality classes termed 
Anomalous Directed Percolation. The analysis of the 
stochastic models reinforce the analogy between epidemic 
spreading mediated by unrestricted Levy fiights and the 
examined STs for a twofold reason: on one side the nu- 
merically evaluated indexes are almost identical to the 
ones found in Ref. j40i] for a lattice model reproducing 
anomalous DP, on the other hand an effective Langevin 
equation has been derived coinciding with that proposed 
for ADP. It is worth stressing that it is highly nontriv- 
ial that the deterministic systems here investigated ex- 
hibit scaling properties in quantitative agreement with 
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those found for stochastic models, hke the lattice model 
studied in [13] and the RM-model. Moreover, the accu- 
racy achieved in the investigation of the STs for the the 
stochastic RM-model suggests that these models can rep- 
resent a valid alternative to the use of lattice dynamics for 
the investigations of non-equilibrium phase transitions. 

The study of smooth continuous maps has revealed a 
behaviour of the critical exponents similar to the one ob- 
served for ADP, namely the indexes vary with continuity 
with the power exponent cr, albeit their values are differ- 
ent from those found for ADP. Moreover, the estimated 
exponents do not correspond to any known universality 
class and this represents a challenge for future theoret- 
ical investigations. In particular, the natural candidate 
to explore is the field equation ([7]) with noise amplitude 
g{n) oc n in order to understand if an Anomalous Multi- 
plicative Noise class could be defined. 
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APPENDIX A: FIELD DESCRIPTION FOR THE 
RM MODEL 

The aim of this appendix is to provide a self-contained 
heuristic derivation of the field equation associated to 
the ST for spatially extended systems with power-law 
decaying interactions, where the transition is controlled 
by finite amplitude effects. In particular, we focus on 
the Random Multiplier model (H]) for A = 0, which 
closely reproduces the ST for discontinuous (or quasi- 
discontinuous) coupled maps. In our derivation we fol- 
low Ref. ,17J where, by introducing a suitable spatio- 
temporal coarse-graining, it has been shown that the 
nearest-neighbor version of the model can be "effectively" 
described by the RFT associated to ordinary DP. In the 
following, we show that the power law coupled version 
of the model can be reduced to Eq. which was pro- 
posed in Ref. [i^l to describe ADP in epidemic spreading 
processes mediated by Levy flights. 

For the fully coupled case, the RM model ^ can be 
rewritten as follows : 

= (1 + V2 + V^)u;,(t) (Al) 

where Vi plays the role of ivi, and is the discretized 
Laplacian operator: 

Vlw^it) = ^w,+i{t) + |u),_i(i) - ew,{t) , (A2) 

while represents a discretized fractional derivative, 
which corresponds to the most relevant term in the small 
momentum expansion of the the following discretized 



convolution sum 

M 

V,w,(i) - -— > — . A3) 

' ^ ' m=2 

Notice that the sum already contains a short distance 
cutoff that should be anyway considered to have a mean- 
ingful definition of fractional derivatives [H, [13] . The 
stochastic variable Wi{t) e [0 : 1] evolves according to 
(|4]), and the positive parameter e represents the am- 
plitude of the spatial coupling. The constant ri'{a) = 

2 Sm=2('™) '^"^ is a normalization factor. Periodic 
boundary conditions are imposed. 
Let us formally rewrite Eq. (|4|) as 

w,{t + 1) = 2av,{t) - a^vlit) + g{v)£,[{t) (A4) 

where the term represents a zero-average (5-correlated 
noise term with unitary variance. In order to recognize 
that the above expression recovers the original model it 
is enough to notice that 

m ^ ^^[Uht) ^ {u^.m , (A5) 

^t, being the dicotomic noise term 

y avi(t), w.p. I - p 

whose average (^^) and variance g'^{v) have the following 
expressions: {^v) — {2a—l)v—a'^v'^ a.nd g^{v)^avSa^v^+ 
Sa^w^-a'^w'* (for details see [TTj). 

We can now introduce a coarse-grained variable 
n(x,t) = Wi{t) (where the bar denotes an average over a 
suitable space-time cell), in terms of which (jAip can be 
written as 

v,{t) = n{x, t) + ^V'^n[x, t) + ^^71(0;, t) 

where the constant c(e) takes in account the presence of 
the cutoff and various normalization factors. The coarse- 
grained evolution equation is then derived from (|A4[) and 
reads 

dtn{x. t) = (2a - \)n{x., t) + aeV'^n{x, t) + ac{e)'V''n{x, t) 
~a^n^{x, t) - —{eV'^n[x, t) + c(e)V''n(x, t)f 
-a^n{x, t){eV'^n{x, t) + c{e)Vn{x, t)) 
+5 (^{x, t) + I V2n(x, t) -f ^V"n(x, t)^ p{x, t) (A7) 

where the coarse-grained noise term p(a;, t) is Gaussian 
and space-time 5 correlated. In proximity of the transi- 
tion the terms of order (V^n)^, (V^n)^, nV'^n, and nV^n 
can be shown to be irrelevant, and also the terms VV^n 
and VV'^ri entering in the noise amplitude g{. . .). 
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By discarding the irrelevant terms one obtains: 

dtn = ((2a- l) + a£V2 + ac(e)V'")n(a;,t) 
—a^n^{x. t) + y^arip{x, t) 



that is essentially the same Langevin equation proposed 
to describe anomalous DP in [33, EOl • 
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